Emission spectrum of quasi-resonant laterally coupled quantum dots 
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We calculate the emission spectrum of neutral and charged excitons in a pair of laterally coupled 
InGaAs quantum dots with nearly degenerate energy levels. As the interdot distance decreases, a 
number of changes take place in the emission spectrum which can be used as indications of molecular 
coupling. These signatures ensue from the stronger tunnel-coupling of trions as compared to that 
of neutral excitons. 

PACS numbers: 73.21. La, 78.67.Hc, 78.55.Cr, 71.35.-y 



I. INTRODUCTION 

There is current interest in developing quan- 
tum dot molecules (QDMs) for quantum information 
processingjii^ Great progress in this direction has been 
made using vertical QDMs fabricated with epitaxial 
growth techniques. These structures are formed by pairs 
of vertically stacked QDs with a thin barrier in between. 
Because the constituent QDs are generally asymmet- 
ric, electric fields applied along the coupling direction 
are used to control the charge and the tunnel coupling 
strength)^""— 

In the last years, developement of QDMs formed by 
laterally coupled QD pairs is also being pursued. ?ri2 In 
these structures, the QD pair is formed inside the same 
epitaxial layers. This architecture is expected to offer 
some advantages as compared to vertical QDMsii^ti^ In 
particular, independent control of QD charging and tun- 
neling strength should be possible through the applica- 
tion of vertical and lateral electric fields, respectively. 
Also, simultaneous coupling of multiple QDMs, which is a 
requirement for the scaling of qubit operations, should be 
feasible through the implementation of individual electric 
gates for each QDM. On the other hand, lateral QDMs 
present their own obstacles. For example, the distance 
between the centers of the QDs is much larger than in 
vertical structures, which implies weaker tunneling. In- 
deed, despite recent advances in the electrical engineer- 
ing of lateral QDMs, which include the demonstration 
of Coulomb blockade chargingi^""— and tuning of optical 
resonancesiS,, the application of lateral electric fields pro- 
ducing clear spectroscopic signatures of molecular cou- 
pling remains a challenge. 

Understanding the differences between the physics of 
lateral QDMs and the better known case of vertical 
QDMs will contribute to determine the conditions for 
optimal device performance, be it vertical, lateral or 
hybridji^ Many of the differences arise from the dis- 
tinct synthetic routes, which lead to characteristic struc- 
tural and compositional profiles i^^'^° Thus, in lateral 
QDMs it has been shown that the weaker quantum con- 
finement of the QDs along the coupling direction leads 
to an enhanced role of the excited orbitals within the 
QDs.— i2^ Also, the nonradiative relaxation dynamics of 
lateral InGaAs QDMs has been shown to follow a singu- 



lar path,^i2^ 

In this paper, we study the optical resonances of neu- 
tral, negatively-charged and positively-charged excitons 
in lateral InGaAs QDMs as a function of the interdot 
distance. The dots have slightly different energies, as is 
often the case with current growth techniquesi^'^ i^^'^^i^^ 
We find that tunnel-coupling of trions is stronger than 
that of neutral excitons. This leads to signatures in the 
photoluminescence spectrum at zero electric field which 
can be used to distinguish uncoupled QD pairs from cou- 
pled QDMs. The convenience of using trions instead of 
neutral excitons for quantum information protocols is dis- 
cussed. 



II. THEORETICAL MODEL 

Our Hamiltonian for excitons and trions can be written 
in the second quantization as: 
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where Ef (E^) is the electron (hole) energy in the single- 
particle state \i) (|p)), ef/ci (hp/hp) is the electron 
(hole) creation/annihilation operator, i.e., {ef e^ . . . 
){h+h+...) \0h)%) = \pq...)\i ]...), and {ij\V\kl), 
{pq\ V \rs) and {ip\ V \qj) are the electron-electron, hole- 
hole and electron-hole Coulomb matrix elements respec- 
tively. Electron-hole exchange term is safely neglected 
because in the absence of magnetic fields and spin mix- 
ing, it has no effect on the observed emission spectrum. 
To calculate the single- particle states, we use a two- 
dimensional effective mass Hamiltonian for electrons and 
heavy holes: 
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where i = e,h is the index denoting electrons or holes, 
m* is the effective mass, Pa — —iKS/a the momentum 
operator, and V{x, y) the quantum confinement poten- 
tial. For a pair of QDs whose centers are separated by a 
distance d, V{x,y) is defined with intersecting parabolic 
potentials: 
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Here ujij^ i'-^i.ii) is the confinement frequency in the left 
(right) QD. A contour plot of this potential profile for 
a slightly asymmetric QDM is shown in Figure [T] Sim- 
ilar models were previously adopted to study laterally 
coupled QDs2i^. 
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FIG. 1. (Color online). Contour plot of the confinement po- 
tential for slightly asymmetric, laterally coupled QDs with an 
interdot distance d=40 nm. 



where V is the polarization operator—, E'^ is the energy 
of the state \s) and F(£') is a Lorentzian curve centered 
at energy E, which simulates the intrinsic bandwidth of 
the transition. Pi(T) is the thermal population distri- 
bution function for the initial state at temperature T. 
We shall consider two possible dynamics for the emis- 
sion process: (i) emission takes place from all the excited 
states, which are populated through non-resonant excita- 
tion, prior to their thermal relaxation. We simulate this 
adopting the same occupation probability for all the ini- 
tial states, Pi(T) = 1; (ii) recombination once the thermal 
equilibrium is reached. In this case we assume a Boltz- 
mann distribution Pi{T) = Z — e~~Er, Here gt {go) is 
the degeneracy factor of the state \i) (ground state), AEi 
the energy difference between \i) and the ground state, 
k the Boltzmann constant and Z the normalization con- 
stant. 

In our calculations we consider laterally coupled 
InGaAs/GaAs QDs. We shall use parameters which, 
for long interdot distances, fit the emission spectrum 
observed in single QDs fabricated upon GaAs nanohole 
templatesjiSii^ Effective masses are ml = 0.0324 and 
TO^ = 0.435, as corresponding to Ino.75Gao.25As, and the 
dielectric constant e = 12. 3j^ The electron confinement 
frequencies are huje,L = 35 meV and fujJe,R = 34 
meV, so that the right QD is slightly bigger than 
the left QD. This is the situation found in some real 
samples.— li^ Hole confinement frequencies are taken 
such that Ih = 0.59 Ze, where li is the characteristic 
length, li — ^jTi/m*uJi. This means that holes are more 
confined than electrons, and their tunneling will be 
negligible. Similar Ih/le ratios have been estimated in 
recent experiments with InGaAs QDs<^ The bandwidth 
of the Lorentzian function in our spectral simulations is 
0.05 meV. For the many-body calculation, we project 
Hamiltonian ([T|) onto all the single-particle configu- 
rations which can be obtained by combining the 24 
lowest-energy electron and hole spin-orbitals. The origin 
of electron (hole) single-particle energies is taken at the 
edge of the conduction (valence) band, disregarding the 
z-confinement energy, which is constant for all the states. 



Hamiltonian ([2|) is integrated numerically using a fi- 
nite differences scheme. The resulting electron and hole 
states are plugged into Hamiltonian ([1]), together with 
the Coulomb matrix elements, which we calculate using 
statistical methods. The Hamiltonian is then solved us- 
ing a configuration interaction routine<^ In this way, we 
obtain an accurate description of the ground and excited 
states of the QDM system. To simulate their emission 
spectrum, we use the dipole approximation and Fermi's 
golden rule;^- The recombination probability from an ini- 
tial few-body state \i) to a final state |/) with one less 
electron- hole pair, at an emission frequency uj is then 
given by: 



^{uj) ex |(/| P |z)|2 T{E^ ~ Ef - hLo) p,{T), (4) 



III. RESULTS AND DISCUSSION 

A. Exciton and negative trion 

Figure [5] shows the emission spectrum of the QDM as a 
function of the interdot distance d, from d = 40 nm (iso- 
lated QDs) to d = 20 nm (strongly coupled QDs). Solid 
red and dashed green lines are used for negative trions 
(AT") and neutral excitons (X"), respectively. Here we 
assume an out-of-equilibrium system. In this way, we 
simulate optically active transitions from ground and ex- 
cited states of both QDs, as normally observed in the 
photoluminescence spectra of verticalA"— and laterali^"— 



QDMs. The most remarkable result of Fig. [2] is that the 
X^ and X^ resonances follow very different evolutions 
as the QDs are brought closer together. To understand 
the different behavior, we need to analyze each excitonic 
complex in detail. 
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FIG. 2. (Color online). Emission spectrum of negative trions 
(solid red lines) and neutral excitons (dashed green lines) for 
different interdot distances. Grey lines are a guide to the eyes. 
See text for the meaning of labels. 

In the limit of the uncoupled QDM, d = 40 nm, X° 
shows two resonances of similar intensity split by '^ 1 
meV. As d is reduced, the two resonances are slowly red- 
shifted, but the energy splitting remains approximately 
constant. This behavior can be understood from the ex- 
citon states, plotted in Fig. [3l The emission resonances 
correspond to the two lowest X'^ states. At d = 40 nm, 
these are the direct excitons, with the electron-hole pair 
in the big (low-energy resonance) and small (high-energy 
resonance) QD, see right-side insets in the figure. Indirect 
excitons (electron and hole in opposite QDs) are much 
higher in energy and optically dark because the electron- 
hole overlap is negligible. As the interdot distance de- 
creases, tunnel coupling allows electrons and holes to 
start hybridizing and form incipient molecular orbitals, 
even if the QDs are not exactly resonant. The lowest 
X'^ state corresponds to bonding (nodeless) electron and 
hole orbitals, so it is redshifted. The first excited state 
corresponds to bonding electron but antibonding hole or- 
bital (note the presence of a node in the green line of the 
upper- left inset). Yet, this makes little difference because 
hole tunneling is negligible. As a result, the two states 
are redshifted by a similar amount. 

The X^ emission spectrum shown in Fig. [2] differs 



from previous calculations for asymmetric lateral QDMs, 
where the two lowest X'^ states were at some point 
blueshifted by Coulomb coupling J This is due to the dif- 
ferent confining potentials in the models. In Ref. 0, the 
interdot barrier was a step-like potential. Instead, we 
use a parabolic potential such that, with decreasing d, 
not only the QDs are closer but also the barrier is lower. 
For this reason, tunnel coupling effects rapidly overcome 
Coulomb coupling in our system. Real QDMs are lens- 
shaped QDs with some degree of lateral contact iiS As the 
dots are brought closer, the height in the contact region 
increases. We then believe our model is more appropri- 
ate. 
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FIG. 3. (Color online). Energy of the lowest exciton states 
as a function of the interdot distance. The insets show the 
charge density of the electron (red line) and hole (green line) 
within the exciton complex for the optically active states, at 
d = 40 nm (right) and d = 20 nm (left). 



Next we investigate the X~ spectrum. As can be seen 
in Fig. [21 at d = 40 nm there are four resonances associ- 
ated with X^ states. Two resonances -(I) and (2)- are 
a few meV below those of X^, while two are degener- 
ate with the X^ ones. As d is reduced, several features 
show up which differ from the simple behavior of neutral 
excitons: (i) the number of resonances increases (up to 
10); (ii) the two low-energy resonances invert their or- 
dering; (iii) the degeneracy between the high-energy X~ 
and the X'^ resonances is lifted. Significantly, most of 
this occurs even before the neutral exciton starts feeling 
tunnel coupling (see e.g. d = 34 nm). 

To understand the X~ emission spectrum, we study 
the energy levels of the initial {X~) and final (le~) 
states, which are plotted in Fig. 21 The electron. 
Fig. mjb), displays a simple dissociation spectrum typ- 
ical of diatomic heteronuclear molecules. At d = 40 nm 
the electron is in one of the two QDs (see right insets). 
With decreasing d, it forms bonding and antibonding 
(noded) molecular orbitals (see left insets). As for the 
trion. Fig. mja), the two lowest-energy states correspond 
to the so-called direct trion states. Their behavior is the 
same as that of X^. Namely, at d = 40 nm all the carri- 



ers are localized inside one of the QDs, and they are red- 
shifted with decreasing d owing to electron hybridization 
(see insets). As indicated by the vertical arrows in Fig.lH 
resonance 1 (2) in Fig. [5] originates in the transition from 
the direct X~ state in the right (left) QD to the elec- 
tron state in the same QD. The energy splitting between 
the two initial states is larger than that between the fi- 
nal states. As a result, transition (1) is lower in energy. 
However, as d decreases the situation reverses. Transi- 
tion (1) involves bonding initial and final states, whose 
energetic stabilizations due to tunnel-coupling compen- 
sate each other. As a result, the energy of the transition 
is little sensitive to the coupling (see Fig.[2|). In contrast, 
transition (2) involves a bonding initial state but an anti- 
bonding final state, so its energy rapidly decreases. This 
is also seen in Fig. [5] and leads to a reversal of resonances 
(1) and (2) at d ^ 32 nni. What is more, because of 
the QDM asymmetry, transitions between bonding and 
antibonding states are not symmetry-forbidden. Then, 
when electron hybridization takes place new transitions 
start developing, (!') and (2'). 

Resonances (3), (4), (5) and (6) in the X^ spectrum 
originate in a set of excited trion states with energy 
~ 0.055 eV at d = 40 nm. These are the so-called in- 
direct trions, where one electron sits in a QD and the 
remaining electron-hole pair sits in the opposite QD (see 
insets in Fig. 01). When the electron- hole pair recombines, 
the electron in the opposite QD remains. At d = 40 nm, 
the Coulomb interaction between the two QDs is so weak 
that the indirect trion resonances have the same energy 
as those of X^. With decreasing d, however, the behavior 
of indirect trions becomes noticeably different from that 
of neutral excitons. The indirect trion states split into 
two groups. One group is formed by states where the 
two electrons form a singlet spin configuration. These 
states increase in energy and eventually anticross with 
higher states. This destabilization is mainly caused by 
the gradual enhancement of the antibonding character. 
In contrast, the group of states with electrons forming a 
triplet remains almost insensitive to tunnel-coupling due 
to the spin-blockade. This effect is analogous to that 
observed in vertical QDMs under electric fields, where 
the Pauli exclusion principle prevents triplet states from 
tunneling.^ The characteristic evolution of indirect tri- 
ons vs. d is responsible for the lifting of the degeneracy 
with X'^ resonances observed in Fig. [21 Last, as in the 
case of direct trions, the QDM asymmetry enables transi- 
tions from indirect trions to either bonding or antibond- 
ing electron states, which gives rise to new resonances 
when electrons hybridize, (5') and (6'). 

The different evolution of X~ and X" resonances in 
the emission spectrum of QDMs can be used to elucidate 
the presence of molecular coupling. The appearance of 
additional trion resonances should be a clear signature, 
but they may be too weak for typical interdot distances 
{d = 30 — 40 nm in Ref . llOl ) . A more evident signature 
should be the splitting between X° and indirect X~ res- 
onances, which involves intense optical resonances and is 
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FIG. 4. (Color online). Energy of the lowest negative trion 
states (a) and electron states (b) as a function of the interdot 
distance. The lateral insets show the charge density of the 
electron (red line) and hole (green line) at d = 40 nm (right) 
and d = 20 nm (left). Vertical arrows indicate the transitions 
observed in Fig. [2l 



visible for long interdot distances {d ^ 36 nm in Fig. [JJ . 
Finally, the reversal of resonances (1) and (2) may also 
be noticed in experiments. For d < 34 nm, this im- 
plies that trion emission from the bigger QD occurs at 
higher energy than that from the smaller QD, contrary 
to the ordering of A"" resonances. As a matter of fact, 
this anomalous spectral ordering has been recently found 
in photoluminescence measurements of lateral QDMs J^ 



B. Positive trion 

We next study the emission spectrum of the positive 
trion (A+). In this case, assuming equal population for 
all the initial states would lead to a complex spectrum 
with many resonances involving highly-excited initial and 
final states. Since we are interested in the main transi- 
tions only, the emission spectrum plotted in Fig. |5l has 
been computed assuming thermal equilibrium at T = 80 
K. At this temperature, the A"+ states with significant 
population are essentially the same as previously studied 
for X~ . The only remaining signature of transitions be- 
tween excited states is the low-energy resonance labeled 
with an asterisk in Fig. [5l which will not be studied here. 

Comparing Figs. [5l and [H one can see that the A+ 
emission pattern differs considerably from the X~ one. 
To analyze the emission spectrum, in Fig.lHlwe represent 
the initial (A+) and final {lh~^) states as a function of 
the interdot distance. The first important difference with 
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appearance of additional resonances, (1') and (3'), at rel- 
atively long interdot distance (see e.g., d — 36 nm). 
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FIG. 5. Emission spectrum of positive trions for different 
interdot distances. Dashed lines are a guide to the eyes. 



respect to X^ can be seen in the final state. Fig. [BJb). 
Unlike for electrons, the two lowest hole states remain al- 
most unaltered throughout the entire range of d, due to 
the negligible tunneling. This means that the evolution 
of the X'^ emission spectrum can be entirely interpreted 
from the initial X+ states. The second important differ- 
ence is in the X+ states, plotted in Fig. [BJa). There we 
find two direct and two indirect states, as in the X~ case, 
but now their energy ordering is reversed. At d = 40 
nm, the two lowest states are indirect trions, with one 
hole in each QD and the electron in either the left or 
right QD (see insets). We note that indirect X+ ground 
states have also been obtained in 3D atomistic calcula- 
tions^^, and are due the larger magnitude of the hole-hole 
Coulomb repulsion as compared to the electron-hole at- 
traction, i.e. \Vhh\ > \Veh\- The situation is the opposite 
to that oi X~, where |Ke| < [K/il, which favors direct 
trion configurations. 

The resonances labeled as (2) and (4) in Fig. [5] orig- 
inate in the recombination of the direct trions. As the 
interdot distance is reduced they redshift slowly. The 
evolution of the indirect trion resonances, (1) and (3), 
is however drastic. The two resonances become rapidly 
split with decreasing d, reflecting the strong splitting of 
indirect trions in Fig. IHKa). Such a large splitting follows 
from the strong tunneling of indirect X^. Unlike for X'^, 
the electron in indirect X'^ can tunnel back and forth 
without losing the strong on-site Coulomb attraction Veh ■ 
This clearly favors molecular coupling. "^^ Furthermore, 
the electron delocalization over the entire QDM enables 
effective recombination with both holes, leading to the 
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FIG. 6. (Color online). Energy of the lowest positive trion 
states (a) and hole states (b) as a function of the interdot 
distance. The lateral insets show the charge density of the 
electron (red line) and the hole (green line) at d = 40 nm 
(right) and d = 20 nm (left). Vertical arrows indicate the 
transitions observed in O 



C. Low^ temperature spectra 

In this section, we re-examine the emission spectra as- 
suming thermal equilibrium at low temperature. This 
greatly simplifies the spectrum, as only populated initial 
states contribute to the spectrum. We consider T = 25 
K, which is enough for both QDs to have finite occupa- 
tion. Fig. [7] shows the emission spectrum for X~ (solid 
red lines), X" (dashed green lines) and X^ (dotted blue 
lines). The evolution of X° and X~ {X~^) resonances is 
essentially that of the direct (indirect) species discussed 
above. Under these conditions, the only signatures of 
coupling remaining are: (i) the appearance of weak X~ 
resonances, labeled in red as (!') and (2'), and (ii) the 
presence of strong low-energy X^ resonances. These 
transitions arise from the ground (indirect) X~^ state, 
which soon split into (1) and (!') and undergo a rapid 
redshift with decreasing d. Thus, at d = 40 nm we start 
with a typical spectrum of isolated QDs, where X° and 
X^ resonances are close in energy, while X~ resonances 
are a few meV below.— i^^ As d is reduced, the X"*" reso- 
nances redshift and approach those oi X~ . 
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FIG. 7. (Color online). Emission spectrum of negative trions 
(solid red lines), neutral excitons (dashed green line) and pos- 
itive trions (dotted blue lines) for different interdot distances. 



interdot distances than X°. The signatures include: (i) 
the appearance of additional optical resonances; (ii) the 
inversion of the ordering of direct trion resonances cor- 
responding to the big and small QDs; (iii) the lifting of 
the degeneracy between indirect trion and neutral exci- 
ton resonances; (iv) a pronounced redshift of the indirect 
X+ resonance, which approaches that of X^ . We note 
that signature (ii) has been actually observed in recent 
photoluminescence experiments^^ We have also shown 
that X~ and X+ develop different emission patterns. 

These results are particularly valuable in the view of 
recent experimental progress to control the charge of ex- 
citonic complexes in QDMs.— "— It has been recently 
shown that trions are more sensitive than excitons to elec- 
tric fields, owing to their net electric charge.— This, to- 
gether with our results, suggests that current attempts to 
utilize lateral QDMs should focus on trion species rather 
than neutral excitons. In both positive and negative tri- 
ons electron tunneling can be efficient. Yet, positive tri- 
ons might be less suited because the spin of holes is very 
sensitive to the in-plane confinement anisotropy."^^'^^ In 
a lateral QDM, the X+ ground state contains one hole 
in each QD. The holes then feel the mutual Coulomb 
repulsion as a source of anisotropic confinement, which 
translates into spin mixing. This could hamper the im- 
plementation of spin-based quantum information proto- 
cols. 



IV. CONCLUSIONS 



We have compared the emission spectrum of X", X~ 
and X+ in asymmetric lateral QDMs as a function of 
the interdot distance. We have shown that tunneling of 
trions is stronger than that of neutral excitons. This is 
because the net Coulomb interactions the electron has 
to overcome in order to tunnel are smaller. As a result, 
X display signatures of molecular coupling at longer 
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